Cooperative Fluorescence from a Strongly Driven Dilute Cloud of Atoms 
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We investigate cooperative fluorescence in a dilute cloud of strongly driven two-level emitters. Start- 
ing from the Heisenberg equations of motion, we compute the first-order scattering corrections to the 
saturation of the excited-state population and to the resonance-fluorescence spectrum, which both 
require going beyond the state-of-the-art linear-optics approach to describe collective phenomena. 
A dipole blockade is observed due to long range dipole-dipole coupling that vanishes at stronger 
driving fields. Furthermore, we compute the inelastic component of the light scattered by a cloud 
of many atoms and find that the Mollow triplet is affected by cooperativity. In a lobe around 
the forward direction, the inelastic Mollow triplet develops a spectral asymmetry, observable under 
experimental conditions. 
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Experimental progress in controlling light-matter in- 
teraction, e.g., in cold atomic clouds and solid state de- 
vices, has in recent years given rise to several propos- 
als and demonstrations of using collections of atoms for 
quantum- information processing 0, When strongly 
driven, a single two-level emitter exhibits a spectral 
triplet, the so-called Mollow triplet which, e.g., has 
been used for generation of heralded single photons and 
entangled photons from solid state quantum dots [4J. In 
collections of many emitters, cooperative phenomena in- 
duced by interatomic dipole-dipole interaction have been 
predicted for weak or no driving leading to cooperative 
decay rates 0, @| and modified Lamb shifts @, [1] that 
have been observed experimentally 0-ll|. While the 
combination of large collections of emitters and strong 
light-matter interaction is surely realizable in the labora- 
tory, theory is faced by a challenge; this setting is difficult 
to tackle theoretically due to the complex nature of the 
non-linear many-body problem. 

State-of-the-art quantum-electrodynamics theory of 
driven atomic clouds typically consider the decay of ini- 
tially inverted systems [a, M , sing le-photon excitations 
in the many-atom case [a, I12h18| , few strongly driven 
atoms 1|J- 21 1 , or interference of light emitted by strongly 
driven non-interacting atoms [22| . 

In this Letter, we report how interatomic interactions 
influence the saturation of the excited-state population 
and the cooperative fluorescence spectrum of a strongly 
driven cloud of two-level atoms. Surprisingly, our results 
show that even when the dipole-dipole interaction be- 
tween any pair of atoms is weak, such as in dilute clouds, 
the collective interatomic coupling is important for the 
nonlinear response. 

Model. — We consider N identical two-level atoms, 
where the m th atom at position R m has a ground state 
\g m ), and an excited state |e m ), separated by the tran- 




FIG. 1. (Color online) (a) Sketch of the setting. A homoge- 
neous laser field of frequency ojo is incident in the direction of 
wavevector ko onto a spherical Gaussian-distributed cloud of 
size cr r containing TV identical atoms. Atoms at positions R m 
and R„ are coupled via the radiation field through G mn • (b) 
Forward-directed elastic part of the fluorescence spectrum di- 
vided by N 2 , calculated using Eq. ((5} with N — 3 (blue dash- 
dotted line), N = 30 (red solid line), and N = 30.000 (green 
dashed line) and fixed bo = 0.1 pumped with A/T = —2.5 
and Q-r/T = 5. The emitted light is detected in the far field 
in the direction of the wavevector of the scattered light ka at 
an angle 9d- 



sition energy hui a - The atoms are driven by a plane- 
wave laser with wavevector ko, frequency ujq, and am- 
plitude Eo, see Fig. HJa). The full Hamiltonian has 
the form H = Ha + Hf + Hi- Here, Ha is the 
free atomic Hamiltonian Ha = ^ ra fiw a ^(t), where 

^m(*) = l(\ e m) (e m \ - \g m ) (5m |) is the population- 
inversion operator of the m th atom. The free-field Hamil- 
tonian is Hf = J2x ^ w Aa^(i)aA(i), where d\(t) is the 
bosonic annihilation operator of the photonic mode A 
with frequency uj\. The electric-dipole Hamiltonian Hi = 
— J2 m f J 'm(t) ' E(R m ,t) describes the light-matter inter- 
action. Here (i m {t) = fj,^ n S+(t)e iuot - ik °- Rm + H.c. is the 
dipole operator, S+(t) = \e m ) (g m \ exp(ik • R m - iuj t) 
the raising operator of the m th atom rotating in the frame 
of the incident field, S^(t) = [S+(t)y the corresponding 
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lowering operator, and fi m = (g m |/x m | e m ) is the dipole 
moment. Finally, E(R m ,t) = i J2x 9\ e \e lkx ' Rm a\(t) + 
H.c. is the electric-field operator with e.\ the polarization 

vector of mode A and g\ = J where eo is the vacuum 
permittivity. 

Dynamics. — We work in the Heisenberg picture and 
after the Born-Markov approximation arrive at the tarna- 
tions of motion for the operators of the m th atom [23] 



dt m 



-(T/2 + iA)S^ 



dt 



qz 



-r s. 



+ i 



2*£ G 

1\ ifln 
+ 2 + — 



G mn S m S n 



q+ qZ 



q+ _ q- 



H.c. 



■F* 



(la) 



(lb) 



Here, SIr = |/x m • Eo|/7i is the Rabi frequency, T = 
4^i 2 o;^/(3?ic 3 ) the spontaneous-decay rate, fi — |/i m | is 
the magnitude of the dipole moment that is equal for all 
atoms, c the speed of light in vacuum, A = luq — w a — rj 
the detuning between the driving field and the atomic 
resonance, and r\ the Lamb shift. Eqs. (Q} are derived 
in the rotating-wave approximation for the atomic oper- 
ators. However the counter-rotating terms in the inter- 
action Hamiltonian are maintained in order to make a 
useful connection with classical optics 14, 2^, 24 1, since 
then the dipole-dipole coupling terms G mn are related to 
the classical Green tensor G [25j by 
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where R mn = R m - R„, i? m „ = |R„ m |, fc a = w /c, and 
we use the scalar model for G that is justified for dipole- 
dipole coupling in the case of dilute clouds Finally, 
the terms F + and F z in Eqs. ([T]) are the Langevin oper- 
ators, which are given by normal-ordered combinations 
of products of atomic and field operators [26l. l27j|. 

Eqs. (JTJ) describe the quantum nonlinear dynamics of a 
cloud of atoms driven by a plane wave of light. For a sin- 
gle atom, these equations reduce to the well-known opti- 
cal Bloch-Langevin equations (26| . Another simple limit 
of Eqs. (p} that does not suffice for the present work is the 
linear-optics limit. The linear dynamics of Refs. 
flil | is obtained from Eqs. ([TJ by the usual approximation 
Sin — — |, valid for weak driving (fiR/T -C 1) that 
maintains the atoms mainly in their ground states. Fur- 
thermore, by letting Or = 0, Eqs. ([TJ also describe the 
dynamics of initially inverted systems leading to super- 
fluorescence as investigated in Refs. Here we focus 
on nonlinear quantum cooperative effects due to strong 
driving. 



Approximate solutions and validity for dilute clouds. — 
Solving Eqs. ((TJ) for the expectation values scales as A N 
so that exact numerical computations for clouds having 
say N ^S> 100 are beyond reach. In the following we fo- 
cus on dilute clouds and aim for accurate rather than ex- 
act dynamics. This allows the simplifying approximation 
that the dipole-dipole coupling between any two atoms 
is small, i.e., G rn n is treated as a perturbation to first 
order. This approximation greatly simplifies the prob- 
lem and allows for analytic expressions for the expec- 
tation values of single-time operators. Furthermore, for 
two-time correlations the Langevin terms contribute neg- 
ligibly and two-time dynamics can thus be reduced by 
the quantum regression theorem to single-time dynam- 
ics. Some details of the method and calculations can be 
found in Ref. 27j. The approximate solutions are valid 
for small optical thickness, bo = 3N/(koa r ) 2 <C 1. This is 
a more severe restriction than on linear theories 
[lH that are valid for small off-resonance optical thick- 
4A 2 /r 2 )<l 



17, 181. This difference 



ness 6a = bo/ (1 
in range of validity can be understood by the fact that, 
contrary to the linear theories, our approach takes all 
frequencies into account and thus there will always be 
some part of the spectrum which is in resonance with 
the atomic transition energy. 

Steady-state population. — Let us first calculate the 
steady-state population of the m th atom, n m = (S m ) + 
1/2. Let n m = + ni\\ where n^m = s/[2(l + s)] is 
the usual single-atom population, expressed in terms of 
the saturation parameter s = r2 R /[2(r 2 /4 + A 2 )], and 

rim is the first-order correction due to the dipole-dipole 
interactions. Solving Eqs. ([I} as a matrix equation to 
first order in G„ m , we arrive at [271 ] 



(Im{G m }r/2 + Re{G m }A)s 
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where G m = J2 n ^mGmn- There is an interesting 
connection between G m , and the cooperative decay 
rate, Tn, and Lamb shift, t]n'. By averaging over 
atomic positions (denoted by an overbar) and consid- 
ering a spherical Gaussian-distributed atomic cloud of 
root-mean-square size cr r (corresponding to atoms in 
a harmonic potential), we obtain for koo~ x 3> 1 that 
Im{G m } = T{N- rj/[2(2/e cr r ) 2 ] and Re{G m } = T(N - 
l)/[2v / 7r(2fc o'r) 3 ] |27j |. These are respectively Ijy and 
r]jy e.g. found from single-photon scattering [y, [l5j and 
t]n also from the scattering correction to the expectation 
value of the Hamiltonian [3, H| . 

For fcoCr 1, i.e., large clouds, Tn/t]n = y/irkoo- T 
such that tjn is negligible and we obtain as a main result 
that the ensemble-averaged mean excited-state popula- 
tion, n = jjY^m n m = ™ (0) 
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6 a s 
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expressed in terms of the off-resonant optical thickness 
6 A ■ Eq. ((4]) shows that for a Gaussian cloud n« is always 
negative, in other words the dipole-dipole interactions 
decrease the steady-state population. This can be inter- 
preted as a cooperative dipole blockade implying that the 
presence of other atoms in the cloud leads to a less effi- 
cient excitation of the emitters. The effect is illustrated in 
Fig. [21 showing that the steady-state excited-state popu- 
lation decreases with increasing optical thicknesses. The 
nonlinear monotonous increase of the population with s 
illustrates that stronger driving makes dipole-dipole in- 
teractions less important relative to the interaction with 
the driving field, and for £7r >Twe recover the steady- 
state population of noninteracting atoms. This agrees 
with and generalizes theoretical observations for two and 
three atoms lil. 21 1. 

To further corroborate our results, we show in the inset 
of Fig. [2] that Eq. (j4j , agrees with existing single-photon 
multiple-scattering theory [3] in the limit of weak scat- 
tering. In more detail, to lowest order in s and bo Eq. (|4l) 
becomes n/n*- ' = 1 — ^p, in agreement with Ref. |18| . 
Thus, Eq. Q unifies both the known dipole-blockade ef- 
fect for weak driving and the novel inclusion of saturation 
effects of the dipole-blockade for strong driving. 

It is interesting to note that, while the large Gaussian 
cloud considered here always results in a blockade effect, 
< 0, an enhanced population due to cooperative 



coupling could be obtained by either of two ways: i) If 
tjnA dominates over TnT/2 a transition from negative to 
positive first-order correction is obtained by varying 
the detuning, ii) If TnT/2 is negative. The case i) is, 
e.g., obtained for a Gaussian cloud when |A|/y / 7rT > 
fcocr which for a cloud-size of fco<7 r ~ ^0 would need a 
detuning of A/r ~ —100. The case ii) could be obtained 
by controlling the atomic positions (e.g., with an optical 
lattice) since the real and imaginary parts of G mn both 
oscillate around zero as a function of interatomic distance 
and thus careful positioning could give a negative Tn = 
Im{G m }. 

Fluorescence spectrum. — Next as our main investiga- 
tion we study the effect of the dipole-dipole interactions 
on the steady-state resonance-fluorescence spectrum of 
the atomic cloud. By assuming that non-scattered light 
is filtered out, we can write the far-field spectrum at de- 
tection angle 8d as 



S(6 d ,uj)/S = 
Re { hm f°° dr (s+(i + r)S~ (t)) e i4 <— iSk - R ™» 
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with So = k\[^ I (YI-k 2 e\r 2 ) where r is the distance from 
the center of the cloud to the detector, buj = uj — luq, 
<5k = kd — ko, and kd is the wave- vector of photons in 
the detection direction. The spectrum consists of two 
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FIG. 2. Normalized population based on Eq. Q as a func- 
tion of saturation parameter, s, for various values of the 
off-resonance optical thicknesses &a- The circles correspond 
to the small saturation parameters used in the inset. In- 
set: Comparison of n/rv ' based on Eq. (j4]) for weak driving 
(s — 2 x 10~ 12 ) with the linearized single-photon multiple- 
scattering theory of Ref. [ifii |. as a function of off-resonant 
optical thickness &a- 



parts. The terms with m — n in Eq. ([5]) concern photons 
emitted from the TV individual atoms. The m ^ n terms 
correspond to interference between photons emitted from 
different atoms. We emphasize that both intensity and 
interference parts have collective features, as the excita- 
tion of each atom is self-consistently obtained by consid- 
ering the drive by the total field, i.e. the incident field 
plus the field scattered by all the other atoms. We eval- 
uate Eq. ([5]) using the quantum-regression theorem and 
split the result into the elastic spectrum S c i as well as 
the inelastic spectrum Si n , which we discuss separately 
below. 

Elastic spectrum. — Based on Eqs. {l} and ([5]), we cal- 
culate the ensemble-averaged angular-emission pattern of 
the elastic spectrum, S c i(8d). It consists of an isotropic 
part, corresponding to the intensity emission, and a 
strongly forwardly directed lobe, due to the interference 
part of the spectrum [2l} ■ Close to the forward direc- 
tion, Sci(^d) scales as / 2 (#d) = exp{— 2[ko<j r sin(6>d/2)] 2 }. 
The function f(9d) is known from weak-scattering the- 
ory and, e-g^ describes interference in Rayleigh-Gans 
scattering [28J. For clouds larger than the wavelength, 
the forward lobe is the dominant contribution to the 
elastic scattering for detection angles smaller than 8 C — 
2/[ln(iV)fcoO'r], i-e., close to the exact forward direction. 
The forwardly directed emission can be seen in Fig. [ljb) 
where S c i(9d)/N 2 is plotted for different N and fixed bo 
for s « 2. The magnitude of the forward emission is 
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given by [27 1 
S ci (Od = 0)/5 = 



ttN 
l + s 



y(0) 



+ (1 - s)n (1 A 5(5w), (6) 



consisting of a non-interacting part and the first-order 
correction. Interestingly, Eq. © shows how the first- 
order correction to the forward-scattering lobe can be 
expressed in terms of the corresponding first-order cor- 
rection to the steady-state population n^ 1 ' of Eq. Q. 
Surprisingly, the elastically scattered intensity, which is 
proportional to the frequency integral of the elastic spec- 
trum, is not proportional to the atomic population as 
is otherwise found in the linear optics |18j and single 
atom [29| limits. This signifies that detection of scattered 
light is not a direct measure of the atomic population. 

For s ps 2 and 6a ~ 0.004 as used in Fig. QJb) the cor- 
rection to the elastic spectrum due to the dipole-dipole 
interaction is on the order of 10 ~ 5 . While this correction 
appears to be small for the elastic spectrum, we will see 
that in the inelastic spectrum the cooperative effects are 
considerable. 

Inelastic spectrum. — We now turn to the inelastic 
component of the spectrum of Eq. ([5]l. Sin, an d study 
how the Mollow triplet is affected by interatomic interac- 
tions. While some limits of the steady-state population 
and elastic spectrum can be investigated in the linear- 
optics regime, the inelastic spectrum is a truly nonlinear 
quantum optical phenomenon that calls for the theory 
reported in this Letter. 

For non-interacting atoms an angle-independent in- 
elastic emission pattern is found, which is simply N times 
the single-atom Mollow triplet. Interestingly, when inter- 
actions are included also the inelastic spectrum becomes 
angle dependent. This is shown in Fig. [3J where the 
inelastic fluorescence spectra of a cloud of N = 25000 
atoms and size ko<J r — 5000 for strong (fiR/T = 5), off- 
resonant (A/T = —2.5) driving are depicted for several 
detection angles. While the number of atoms N — 25000 
is by far too large for usual numerical calculation meth- 
ods, it is relevant for experimental settings. The cloud 
size corresponds to the experimental value of Ref. [lij ]. 
where a detuning ranging from — 1.9r to — 4.2r was used 
consistent with our A = — 2.5r. 

All spectra shown in Fig. [3] exhibit the typical three- 
peak structure of the single-atom Mollow spectrum, but 
also a spectral asymmetry, which is strongest in the for- 
ward direction, see inset. In contrast, the single-atom 
Mollow spectrum is symmetric, even for off-resonant 
driving. The observed asymmetry depends on the laser- 
atom detuning and gives an increase of the sideband peak 
closest to the bare atomic transition frequency u a . The 
~ 35% enhancement of the peak at the bare atomic fre- 
quency, as shown in Fig. [3j is a result of cooperative 
effects showing the importance of including the dipole- 
dipole interactions when dealing with the fluorescence 
spectrum. 
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FIG. 3. (Color online) Inelastic fluorescence spectrum versus 
rescaled detection angle, fcoo^d, with N = 25000, koa r = 
5000, fiii/r = 5, and A/r = —2.5. Inset: comparison of 
the asymmetric inelastic fluorescence spectrum for 8d — 
(black solid line) with the symmetric non-interacting atoms 
spectrum (black dashed line with green shaded area). The 
vertical red line shows cj a . 



We can now appreciate the pronounced cooperative 
features in the inelastic spectrum as compared to the 
elastic spectrum: for non-interacting atoms the elastic 
spectrum already shows an iV-times enhanced forward- 
directed peak, but the inelastic spectrum does not. While 
the elastic peak is slightly modified due to interatomic 
interactions for the inelastic spectrum the same interac- 
tions create a forwardly-directed peak. In more detail, 
the dipole-dipole interactions create interatomic correla- 
tions such that the forward-directed interference pattern 
is built up in the inelastic spectrum. 

Conclusions, discussion, outlook. — In conclusion, we 
have shown that dipole-dipole interactions in clouds of 
cold atoms affect their optical properties in the strong- 
driving regime, even for dilute clouds. We found analyti- 
cal corrections to the steady-state population and to the 
fluorescence spectrum under strong driving. The anal- 
ysis allows connecting the cooperative decay rate and 
Lamb shift with the Green function governing photon 
propagation. We found that, while a spherical Gaussian 
distributed cloud exhibits decreased atomic excitation, 
also a cooperatively increased atomic excitation is pos- 
sible. Moreover, we have shown that cooperative scat- 
tering persist in the Mollow triplet, which is a hallmark 
of non-classical scattering of light by two-level systems. 
The cooperative effect gives rise to an angle-dependent 
spectrum and is most pronounced in the forward direc- 
tion where it manifests itself as an enhancement of the 
sideband nearest to the atomic-transition frequency. 

While we have considered the simplest model for the 
atoms, the scalar two-level model, it is worth noting that 
the approach used in this work can be generalized, e.g., 
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to account for the full vectorial nature of the atom-light 
scattering and the near-field components of the dipole- 
dipole coupling. We are confident that our results for 
strongly driven dilute clouds will stimulate the study 
of denser clouds where interatomic interactions are ex- 
pected to be even more important. 
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